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Abstract — Digital breast tomosynthesis (DBT) is an emerging 
modality for breast imaging. A typical tomosynthesis image is 
reconstructed from projection data acquired at a limited number 
of views over a limited angular range. In general, the quantitative 
accuracy of the image can be significantly compromised by 
severe artifacts and non-isotropic resolution resulting from the 
incomplete data. Nevertheless, it has been demonstrated that DBT 
may yield useful information for detection/classification tasks and 
thus is considered a promising breast imaging modality currently 
undergoing pre-clinical evaluation trials. The purpose of this 
work is to conduct a preliminary, but systematic, investigation 
and evaluation of the properties of reconstruction algorithms 
that have been proposed for DBT. We use a breast phantom 
designed for DBT evaluation to generate analytic projection data 
for a typical DBT configuration, which is currently undergoing 
pre-clinical evaluation. The reconstruction algorithms under 
comparison include (i) filtered backprojection (FBP), (ii) expec- 
tation maximization (EM), and (iii) TV-minimization algorithms. 
Results of our studies indicate that FBP reconstructed images are 
generally noisier and demonstrate lower in-depth resolution than 
those obtained through iterative reconstruction and that the TV- 
minimization reconstruction yield images with reduced artifacts 
as compared to that obtained with other algorithms under study. 

I. Introduction 

The breast cancer death rate declined by an average of 2.3% 
per year between 1990 and 2002 [1]. This has been attributed 
to earlier detection of breast cancer through mammography 
screening. In a mammographic exam, two x-ray projections 
of each breast are acquired in cranio-caudal and mediolateral- 
oblique views. While mammography is considered the gold 
standard in breast imaging, mass lesion detection is limited by 
anatomic background structures [2], which can obscure lesions 
or mimic a lesion appearance. 

This limitation can be significantly reduced in tomosynthe- 
sis imaging, in which a volume image is reconstructed from 
projection data acquired at a limited number of views over a 
limited angular range. Typically, the number of views ranges 
between 11 and 21, whereas the angle ranges between 15 and 
50 degrees. Tissue structures that are overlaying in conven- 
tional mammography can be resolved in the reconstructed vol- 
ume, and lesions can thus become more conspicuous. Because 
of the data incompleteness in DBT the quantitative accuracy of 
DBT images is generally compromised significantly by severe 
artifacts and non-isotropic resolution. However, the purpose 
of DBT is not to provide an accurate attenuation map of the 
breast. Instead, DBT is aiming only at providing a clinically 



useful image, i.e., to "remove visual clutter" [3], in terms of 
detection/classification tasks. 

The advent of digital detectors for mammography has made 
DBT become tangible since the late 1990s. The development 
of DBT was initialized by a breast mammographer, Daniel 
Kopans and colleagues, at the Massachusetts General Hospital 
[4]. Currently, a number of DBT prototypes have been built, 
implementing a range of acquisition geometries [5], [6], [7], 
[8]. In general, reconstruction algorithms for DBT that have 
been explored fall into two broad categories, namely filtered 
backprojection [5], [6], and iterative methods [8], [9]. 

Accurate image reconstruction in DBT is challenging be- 
cause of the high degree of data incompleteness. Therefore no 
analytic solution exists for quantitatively accurate image re- 
construction from projection data. On the other hand, iterative 
methods have been used to address image reconstruction in 
DBT, and they appear to yield visually improved images over 
those obtained with the non-iterative algorithms. In DBT, high 
in-plane resolution is required in order to visualize microcalci- 
fication clusters. The in-depth resolution of the images is poor 
largely because of the limited angular range. Typically, the 
image sampling is non-isotropic, with an in-plane resolution 
of up to a factor of 10 higher than the in-depth resolution, 
which can lead to artifacts during reconstruction. While this 
anisotropy could in principle be remedied by increasing the 
sampling grid resolution, this approach is not practically 
feasible because of the large data size. Currently, a typical 
tomosynthesis image that adequately covers the breast volume 
consists of about 1.32 s voxels (corresponding to a resolution 
of 0.1 x 0.1 x 1 mm 3 ). 

The purpose of this work is to compare reconstruction 
methods that are currently used for and/or have been proposed 
for DBT reconstruction. Our comparison study is based on 
an analytic phantom that we have developed for the purpose 
of DBT algorithm evaluation. Dominant features of clinical 
relevance have been incorporated into this model, namely the 
overall breast shapes which causes a drop-off in intensity 
towards the skinline, spherical structures to represent tumor 
lesions, as well as microcalcification clusters. 

II. Breast Phantom 

A breast phantom was designed to provide a basis for 
algorithm comparison. It consists of objects with different 
geometric shapes to allow for analytic computation of the 





Fig. 1. The breast phantom used as the basis for comparison of recon- 
struction algorithms in DBT. Images of the breast phantom within 2D planes 
specified by (a) z = 30 mm, by (b) x = 5.5 mm, (c) x = 24.5 mm, and (d) 
x = 41.9 mm, respectively. 



projection data. The overall shape of the compressed breast 
of 50 mm thickness is modeled by a truncated ellipsoid. 
Attenuation coefficients at a photon energy of 30 keV [10] 
were used for phantom materials simulating the fibroglandular, 
tumour and muscle tissues, as well as microcalcifications 
in the breast. In Fig. [T] we display the phantom images 
within 2D planes specified by (a) z = 30 mm, by (b) 
x — 5.5 mm, (c) x — 24.5 mm, and (d) x — 41.9 mm, 
respectively. As shown in Fig. [ij, the crescent- shaped region, 
representing the fibroglandular tissue region of the breast, is 
attached to a rectangular slab of higher attenuation, simulating 
the pectoralis muscle. It can be observed in Fig. [T] that 
numerous test objects of different sizes and contrast levels 
are embedded in the phantom for simulating mass lesions and 
microcalcifications. In particular, a row of three ellipsoids is 
embedded in the pectoralis muscle (see Fig. [TJb)), with equal 
in-plane diameters but varying flatness, thus allowing one to 
evaluate whether shapes of equal in-plane profile, but different 
in-depth profiles, can be resolved. Furthermore, six spheres of 
different diameters, ranging from 5 mm to 15 mm, are equally 
distributed in the fatty area of the breast (see Fig. [ijc)). It 
can also be observed in Fig. [TJd) that three pairs of stacked 
spheres of identical diameter d = 10 mm are embedded in 
the breast region with different spacings of d, 2d, and 3d, 
respectively, for these pairs. Finally, two clusters of small 
spheres of diameters 0.3 mm and 0.15 mm are included to 
model microcalcifications. 

III. Scanning Configuration 

In this work, we consider a tomosynthesis scanning con- 
figuration with parameters similar to those of the first GE 
tomosynthesis prototype for breast imaging [9]. In this con- 
figuration, the source is acquiring data at 1 1 projection views 
uniformly distributed over an arc of 90 degrees. The source-to- 
detector distance is 660 mm, whereas the distance between the 
source and center-of-rotation is 460 mm. The detector plane is 
perpendicular to the line connecting the x-ray source and the 



detector center, at the center (6th) projection, and the detector 
remains stationary while the x-ray source rotates along the 
arc. The detector size is 180 mmx85 mm with the shorter 
side perpendicular to the plane of rotation, and the detector- 
bin size is 200 microns. In our configuration, the object center 
is located 50 mm above the detector. 

IV. Reconstruction Algorithms 

Noiseless projection data were computed analytically from 
the breast phantom for the scanning configuration described 
above. Noisy data were created by adding Poisson noise to 
the noiseless projection data. The photon counts in projection 
data were weighted to account for variations in the source-to- 
detector-bin distance as well as variations in the effective area 
of the detector-bin surface due to oblique x-ray incidence. We 
have used a typical clinical exposure to determine the added 
noise level. From the projection data, images within 3D arrays 
with a voxel size of 0.2 x 0.2 x 1mm 3 were reconstructed by 
use of different reconstruction algorithms. 

A. The FBP algorithm 

For the filtering and backprojection reconstruction (FBP) 
[11], each projection image is filtered with a ramp filter and 
then back-projected using a cone-beam geometry. A von Hann 
window is applied to supress high frequencies. 

B. The EM algorithm 

In each iteration of the expectation-maximization algorihtm, 
the current image estimate fk+i is updated by 



fk+i =fk-M+{-Lr 



} 



(1) 



where g is the projection data, M and M + are the forward 
and backprojection operators, s = ^M, and -^U- indicates a 

M Jk 

division between the vector elements. In our implementation, 
a matched projector/backprojector pair is used with Gaussian 
smoothing between iterations. 

C. The TV-minimization algorithm 

In the TV-minimization algorithm [12], the image is ob- 
tained by solving an optimization problem: 



/* = argmin||/|| Ty such that Mf - g 



(2) 



where / and g are discrete image and data, and M is the linear 
operator representing the cone-beam forward projection and 
/* is the reconstructed image. The parameter e can be selected 
for controlling the impact level of potential data inconsistency 
on the image reconstruction. 

V. Results 

Using the tomosynthesis data produced from the breast 
phantom with the scanning configuration described, we have 
carried out a preliminary numerical studies to compare the 
performance of the FBP, EM, and TV algorithms. In Fig. [2] 
we displayed the images within the plane z = 30 mm 
reconstructed by use of the FBP (left), EM (middle), and TV 




Fig. 2. Images reconstructed within a 2D slice specified by z = 30 mm , 
reconstructed from noiseless data (top row) and from noisy data (bottom row) 
by use of the FBP (left column) EM (middle column), and TV (right column) 
algorithms, respectively. 



(right) algorithms. The true image within the corresponding 
slice is shown in Fig. [ija). As expected, because of the 
strong data incompleteness in tomosynthesis, all reconstruc- 
tions exhibit obvious artifacts. Out-of-plane objects create 
conspicuous ghosting artifacts in all reconstructions, but most 
severe in the FBP reconstruction. We refer to these artifacts 
as structure noise. Among these algorithms, TV algorithm 
appears to produce more uniform images than do others. 
Furthermore, in the FBP reconstruction, edges are "enhanced" 
by the ramp filtering. The noisy reconstructions are visually 
similar, suggesting that the image noise is dominated by 
structure noise caused by out-of-plane objects. It is also ob- 
servable that the noise properties of each algorithm is different. 
The EM reconstruction appears most smooth, while the TV 
reconstruction exhibits some "speckle noise", i.e., spikes, in 
an otherwise smooth image. 

In Figs. [3j and|4| we show the image profiles on three dif- 
ferent lines through the reconstructed image within the center 
plane (i.e., z = 30 mm) of the image. The reconstruction 
properties discussed above can be observed in these profile 
results. Clearly, the FBP reconstruction appears most noisy, 
which is expected because the ramp filtering tends to amplify 
noise. There is a considerable DC shift for FBP reconstruction, 
as well as for the EM reconstruction. From the profile across 
spheres of varying diameters in Fig. [3] it can be observed 




100 200 300 400 500 600 700 800 900 
y [0.2 mm] 



Fig. 3. Profiles on the line specified by x = 24.8 mm and z = 30 mm in 
images displayed in Fig. [2] obtained with the FBP (blue), EM (red), and TV 
(green) algorithms, respectively. Solid curves are for noiseless images, dotted 
curves for noisy images. For comparison, the true profile is depicted by the 
black solid curve. 
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Fig. 4. Profiles on the line specified by x = 34.4 mm. Reconstruction 
algorithms and line types are identical to those in Fig. [3] 



that that in-plane resolution is maintained in DBT, which is 
the reason for the clinical usefulness of DBT. The profile 
in Fig. [4] is on a line through a uniform section of the 
phantom. Uniformity appears to be best reproduced by the 
TV reconstruction. In the TV reconstruction, a shift between 
noiseless and noisy data can be observed. 

We also show in Fig. [5] the images from noisy data within 
planes that are parallel to the source-motion plane and per- 
pendicular to the detector. The true images corresponding to 
those in rows 1, 2, and 3 of Fig. [5] are shown in in Figs.[TJb), 
[TJc), and []Jd), respectively. As the results indicate, none of 
the algorithms recover the true shape of the the ellipsoids of 
varying flatness in the first row. Rather, the variation in flatness 
results in a variation of contrast. Furthermore, an increase in 
overall intensity across the uniform pectoralis muscle can be 
observed for the iterative algorithms. For all reconstructions, 
as shown in the second row of Fig. [5] the spheres of different 
sizes are elongated due to the limited angular range. Objects 
are most uniform in the TV reconstruction and most blurred in 




Fig. 5. 2D slices within images reconstructed by use of the FBP (left 
column), the EM (middle column), and the TV-minization (right column) 
algorithms, respectively. The corresponding true images are displayed in 
Figs. [2]3,[^ and |2jl. 




Fig. 6. Images of a set of small microcalcifications (0.15 mm in diameter) 
within the 2D in-focus plane reconstructed by use of the FBP (left column), 
EM (middle column), and, TV-minimization (right column), respectively, from 
noiseless (row 1) and noisy (row 2) data. 



the FBP slice. The improved in-depth resolution of the iterative 
algorithms is also demonstrated for stacked spheres shown 
in the third row of Fig. [5] Both algorithms clearly separate 
spheres spaced by 2d, while FBP still shows some overlap. In 
all slices obtained through FBP reconstruction, quantum noise 
is quite pronounced. For the iterative algorithms, structure 
noise is dominant. 

Finally, we show in Fig. [6] a region-of-interest (ROI) in 
the in-focus plane of a microcalcification (MC) cluster. The 
top row shows reconstructions from noiseless data, the bottom 
row shows reconstructions from noisy data. The entire cluster 
consists of 6 MC, with another triangular arrangement 6 mm 
below the one depicted. In the noiseless images, patterns 
caused by streaking artifacts can be perceived for the FBP and 
EM reconstructions. Image quality in the TV reconstruction is 
clearly superior in this situation. 



VI. Conclusion 

We presented a systematic comparison of three recon- 
struction methods (FBP, EM and TV), for noiseless and 
noisy phantom data acquired for a tomosynthesis scanning 
configuration. The comparison was based on a phantom that 
was specifically designed for evaluating tomosynthesis images. 
Our study demonstrates that all reconstructed images exhibit 
artifacts, because of the incomplete projection data available 
for reconstruction. However, in-depth resolution was higher in 
images from iterative reconstructions, and therefore ghosting 
artifacts were also less conspicuous in the iterative reconstruc- 
tions. Furthermore, structure noise from out-of-plane objects 
was dominant for iterative reconstructions, while quantum 
noise became visible in the FBP reconstruction. The FBP 
reconstruction also showed artifacts caused by the ramp filter- 
ing. In summary, our study indicates that iterative algorithms 
may provide higher image quality for tomosynthesis image 
reconstruction than FBP-based reconstruction. 
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